Packing of elastic wires in spherical cavities 
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We investigate the morphologies and maximum packing density of thin wires packed into spherical 
cavities. Using simulations and experiments, we find that ordered as well as disordered structures 
emerge, depending on the amount of internal torsion. We find that the highest packing densities are 
achieved in low torsion packings for large systems, but in high torsion packings for small systems. 
An analysis of both situations is given in terms of energetics and comparison is made to analytical 
models of DNA packing in viral capsids. 



PACS numbers: 87.10.Pq, 46.70.Hg, 89.75.Da 

Thin objects are ubiquitous in nature and technology. 
Driven by forces and external constraints, they can un- 
dergo surprisingly complex spatial rearrangements, ob- 
served, for example, in the folding of insect wings in co- 
coons, crumpled wires and paper, or growing tissue [THE]. 
The packing of long, slender objects in cavities emerges 
in many situations in biology and mechanics. It occurs 
in paper jams, when chromatin is stored in the cell nu- 
cleus, DNA is injected into viral capsids [51-fTT]. or when 
endovascular coils are formed in aneurysm surgery [T^ . 
Often, such systems are geared towards high packing den- 
sities. For instance, the amount of DNA packed in a cap- 
sid limits the genetic information the virus can spread, 
and in endovascular surgery, high densities improve the 
long-term stability of the treatment |T3]. 

In this Letter, we consider the packing of a thin wire 
in a rigid spherical container. First, the emergence of or- 
dered as well as disordered packings depending on the 
amount of internal torsion is shown. We show how 
the maximum packing density depends on morphology 
and system size. We find that ordered packings provide 
higher maximum packing densities at large system sizes, 
while disordered ones offer higher maximum densities at 
small system sizes. 

We consider two setups in simulations and experi- 
ments: In the first, straight wires are used that can ax- 
ially rotate at the injection point, thereby having mini- 
mal torsion during the entire packing process. We call 
this the low torsion setup. In the second setup, we aim 
at packings that build up internal torsion. To hinder the 
release of torsion via the free end, we inject precurved 
wires [Hj. While this modification changes the bending 
energy, it surprisingly had no influence on morphology 
and packing density, as torsion is released via the cav- 
ity opening. We thus also suppress axial rotation at the 
opening to arrive at a high torsion setup. 

The experiment consists of a transparent, hollow rigid 
sphere of inner radius R with a small hole to which a 



nozzle is attached. We insert a Nylon or Silicon wire of 
constant cross section radius a by two counter-rotating 
rollers, which allow for a controlled insertion speed and 
large forces. The Young's modulus Y of the wires was de- 
termined from axial extension tests and deflection mea- 
surements ^53^. For the high torsion setup, we use wires 
of constant intrinsic radius of curvature Ri ~ 2R. In 
the low torsion setup, we straightened the wires and al- 
lowed axial rotation between the nozzle and the sphere. 
The packing process is recorded by a camera, and the 
insertion force is measured continuously by a load cell. 
The process stops when the force becomes so high that 
no further insertion is possible. Tomography images are 
taken from the final packing to reconstruct the internal 
structure. 

For the numerical model, we describe the wire by 
its centreline 7, and an orthonormal director field dj, 
i = 1,2,3 that specifies the orientation of the cross sec- 
tions along 7 [TS]. We parametrize 7 by s, < s < L, 
with L the length of the wire. The position of a point on 
the centreline is denoted as r(s). The director field is ori- 
ented such that ds coincides with the centreline tangent, 
r'(s), where ' is short-hand for dg- It follows from the 
theory of space-curves that a vector k exists such that 
d^ = k X d^. k is known as the Darhoux vector and we 
denote its components in the director frame by hi. To 
distinguish between the deformed and undeformed wire 
configuration, we use the superscript °, i. e. the refer- 
ence configuration is described by the smooth curve 7*^, 
directors d^ and Darboux vector k''. If the wire is thin 
compared to its curvature, the internal elastic energy is 
given by 

E^i = \j\s[YA{r',~r','f 

+ YI [(fci - klf + (fc2 - klf] +G{h - kl)^} . (1) 

The first term accounts for axial stretch- 
ing/compression, with the cross section area A = 7ra^. 
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FIG. 1: Morphologies of wires packed into spherical cavities. Top row: The low torsion setup results in an ordered morphology, 
characterized by ring-like coiling (a). A cut through the packing (c) reveils the shell- like inner structure. X-ray tomography 
scans from experimental realizations are shown in (b, d). Bottom row: The high torsion setup produces disordered structures 
(e, g), with corresponding experiments in (f,h). Color represents the bending energy Eb- System parameters are: (j> = 0.23, 
a/R — 0.02 (simulations); <j} = 0.14, a/R — 0.025 (low torsion, experiment); </!> — 0.18, a/R = 0.017 (high torsion, experiment) 



The second term is the bending energy Eb, with 
/ = 7ra^/4 the 2nd moment of inertia. The components 
/cj and ^2 describe the intrinsic curvature of the reference 
configuration. The last term accounts for torsion, with 
G = 2{i+v) being the torsional shear modulus, and 
1/ the Poisson ratio. For naturally straight rods, the 
bending energy simplifies to the well-known expression 
I YlK^ds, where k is the curvature of the centreline. 

For an efficient simulation including self-contact, we 
discretize the wire into N mass-points r^, connected by 
straight edges e^. Tensile energy is modeled by linear 
springs connecting mass-points. For bending and torsion, 
we use the quaternion group to represent the rotation of 
the director frame of each edge, with resulting moments 
and forces given by the gradient of the discretized form of 
Eq.[l][TB]. Contact of the wire with itself or with the cav- 
ity is modeled by a linear repelling force law. We use the 
same stick-slip friction model as in Ref. [2] with friction 
coefficients /is = 0.2 for static and Hd = 0.18 for dynamic 
Coulomb friction - a choice which yields best agreement 
with experiments. No friction is applied to wire segments 
in the injection nozzle for both setups. To prohibit ax- 
ial rotation in the high torsion setup, any resulting axial 
moment in the nozzle is canceled by an equal but op- 
posing one at every time step. Forces and moments are 
integrated in time using a standard predictor-corrector 
method of 6th order, with small viscous damping added 
for equilibration and numerical stability. We fixed the 
wire radius to a = 1 and chose 7 different sphere radii 
R = (4, 5, 10, 13, 20, 40, 50). The Young's modulus 
was set to Y = 5 and ky = 0.001 for the high torsion 



setup, corresponding to an intrinsic radius of curvature 
of Ri = 100. The insertion speed was sufficiently small 
to ensure being in the quasi-static regime, and the sim- 
ulation was stopped when the measured injection force 
reached a fixed threshold. 

The packing process starts with one end of the wire 
being inserted into the cavity. To break symmetry of 
the straight wires in simulations, small random displace- 
ments are initially imposed on all nodes. When the wire 
contacts the cavity walls for the first time, a loop forms 
along a cavity wall. The orientation of this loop is ran- 
dom for naturally straight wires depending on initial con- 
ditions, whereas for curved wires it is in the preferred, 
curved direction. As more wire is inserted, distinct mor- 
phologies emerge for the two setups: 

Low torsion setup: Here, the wire continues to form 
loops which align in parallel with each other. This pro- 
cess leads to ring-like structures with decreasing coiling 
radius. Whenever the coiling radius becomes too small, a 
new coil is started at a different orientation (Fig.[Tj a,b). 
The result is an ordered packing in layers from outside in- 
wards (Fig.[l] c,d) at earlier stages of the process, similar 
to the coaxial coiling model for DNA inside spherical cap- 
sids [2T]. At higher densities, coils prefer alignment per- 
pendicular to the feeding axis, sharing similarities with 
the inverse spool DNA model |11) . 

High torsion setup: In this setup, torsion can only be 
minimized on the expense of bending deformations, lead- 
ing to frequent reorientations of loops. In certain cases, 
torsion becomes so large that figure-eight patterns ap- 
pear, or loops form with smaller radius than imposed by 
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FIG. 2: Top: Maximum packing density 4>max as a func- 
tion of the effective system size a/R for the low torsion (a) 
and the high torsion setup (b). Experimental data are shown 
with empty, numerical data with filled symbols. DNA packing 
densities (V) are added in the low torsion setup for reference 
|17| . Bottom: Length-radius scaling in the high torsion setup 
(c) , and a comparison between the total accumulated internal 
torsional energy Etors as function of (d). 



geometric conditions. As a consequence, the packing is 
disordered (Fig.jlj e,f) and fills the cavity homogeneously 
(Fig. 1 g,h) with a large amount of accumulated torsion 
(Fig.g d). 

We turn to the study of the maximum packing density 
4>max and its dependence on the effective system size a/R. 
4>max follows from the inserted wire length L as 4'max — 
Ltto^ / [A/ZttB?) and is shown in Fig. [2] (a, b) as function 
of a/R for both setups. The system sizes in experiments 
were in the range 0.0069 < a/R < 0.23, with values of 
R[mm] = (6.8, 14.0, 19.9, 22.9, 28.5, 30.3, 38.5, 48.5) 
and a[mm] = (0.16, 0.25, 0.35, 0.4, 0.5, 1.5, 1.93, 2.0, 
3.25). In simulations a/i? was within 0.02 and 0.25. Since 
data points in each setup collapse on a common curve, 
we conclude that (j)max in either setup only depends on 
a/R and not on a and R individually. Furthermore, no 
dependence on Young's modulus Y and on the amount 
of intrinsic curvature was found within the tested range 
R<Ri< 5R. 

In comparison, the low torsion setup (Fig. [2] a) yields 
higher maximum packing densities for large systems 
(small a/R), comparable with experimental data from 
DNA coiling (triangles). As a/R is increased, 4'max in- 
creases slower than (pmax for the high torsion setup. Con- 
sequently, around a/R — 0.2, the high torsion setup 
starts to provide larger (j)max, with packing fractions up 
to 58%. This could be due to the fact that torsion sup- 
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FIG. 3: Rescaled bending energy and injection force (inset) 
of the wire for the low torsion setup (simulations: continuous 
line, experiments: □) and the high torsion setup (simulations: 
dotted line, experiments: o), with a/R= 1/40. The dashed 
line is the prediction of Ref. [llj . Simulation data are aver- 
aged over 9 runs. 



ports the bending of the wire and thus reduces the re- 
quired buckling force, similar to the formation of a DNA 
plectoneme [TS]. 

The strong dependence of 4>rnax on a/R in the high 
torsion setup is captured well by the power-law ^max ^ 
(a/R)", with a = 0.38 ±0.04 (Fig. [2) b). We can further 
conclude that if (j)max scales as a power-law in a/R, then 
the packed wire mass scales as M ^ a^L ^ R^{a/R)°', 
i. e. M/a"" - i^^-a ^^ig. [2] c). The exponent 3 - a = 
2.62±0.04 is close to the experimental value of 2.75 found 
by Gomes et al [19j for the forced crumpling of wires in 
three dimensions. 

The dominant energy in both setups is the bending 
energy E},. To compare results of different cavity radii. 
Fig. [3] shows the dimensionless quantity EhR/YI, which 
yields good agreement between simulations and experi- 
ments !24, . We compare these findings to the packaging 
of DNA in viral capsids. While this process is also ac- 
companied by, e.g., entropic and electronic contributions 
P2i , comparison can be made to the bending-dominated 
DNA model of Purohit et al. [TT]. It predicts 
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in which k = [3dg/(47r^a'*)]^/'^. The only remaining 
parameter is the average segment distance dg, which 
we determined from cross sections of the simulations as 
ds — 2.81a. The prediction agrees well with our data, c.f. 
Fig. |3j dashed line. The dimensionless insertion force 
FR^ /YI obtained analytically from Eq. [2] is, however, 
somewhat smaller than the measurements (Fig. [3] inset), 
due to the fact that we include contact friction. 

The distribution of the bending energy measured at 
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FIG. 4: Bending energy distribution of the packed wire in 
log-linear scale and log-log scale (inset) for a/R = 1/40. Data 
for the high torsion setup (*) can be approximated by a log- 
normal distribution with parameters /i = —0.903 and a = 
1.275 (continuous line). For the low torsion setup (-I-, dotted 
line is shown as a guide to the eye) small bending energies 
dominate due to the ordered coiling. 



4'max gives further insights into the statistical properties 
of the morphologies (see Fig. [4]) . The low torsion setup 
(-I-) features a strong concentration of energies around 
Eh/{Eb) = 0.1 due to the presence of ordered coils with 
large radii. Subsequent coiling rings are squeezed into 
existing ones. They take on a shape reminding of stadium 
racepaths, with almost straight parts that lead to a high 
probability for small Ei,. In contrast, the high torsion 
setup (*) yields a distribution fitted well by a log-normal 
function of the type 
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with /i — —0.903 and a — 1.275. Log-normal distribu- 
tions are usually associated with hierarchical events and 
are found in similar systems of densely packed objects, 
for instance in the ridge-length distribution of crumpled 
paper or two-dimensional crumpled wires [U [51 [5D] . 

To summarize, we investigated the packing of elastic 
wires into spherical cavities using a high torsion and a 
low torsion setup. Low torsion led to ordered packings, 
while high torsion led to disordered structures. We found 
highest packing densities in the first case for large cav- 
ities and in the second case for small cavities, with a 
cross-over as the system size is varied. Our work elu- 
cidates the importance of torsion in the dense packing 
regime, and provides novel insights into the role of the 
system size with relevance for, e.g., the surgical treat- 
ment of aneurysms. The presented results are largely 
independent on Young's modulus, the amount of intrin- 
sic curvature within the tested range, and friction. The 
role of material nonlinearities, however, remains an im- 
portant open question in view of such applications. 
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